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ABSTRACT 

Variability is a property shared by practically all AGN. This makes variability selection a possible 
technique for identifying AGN. Given that variability selection makes no prior assumption about spectral 
properties, it is a powerful technique for detecting both low-luminosity AGN in which the host galaxy 
emission is dominating and AGN with unusual spectral properties. In this paper, we will discuss and 
test different statistical methods for the detection of variability in sparsely sampled data that allow full 
control over the false positive rates. We will apply these methods to the GOODS North and South fields 
and present a catalog of variable sources in the z band in both GOODS fields. Out of 11931 objects 
checked, we find 155 variable sources at a significance level of 99.9%, corresponding to about 1.3% of all 
objects. After rejection of stars and supernovae, 139 variability selected AGN remain. Their magnitudes 
reach down as faint as 25.5 mag in z- Spectroscopic redshifts are available for 22 of the variability 
selected AGN, ranging from 0.046 to 3.7. The absolute magnitudes in the rest-frame z-band range from 
— 18 to -24, reaching substantially fainter than the typical luminosities probed by traditional X-ray and 
spectroscopic AGN selection in these fields. Therefore, this is a powerful technique for future exploration 
of the evolution of the faint end of the AGN luminosity function up to high redshifts. 

Subject headings: astronomical data bases; catalogs; galaxies: active; quasars: general 



1. Introduction 

In the last 20 years or so, findings of strong corre- 
lations between the black hole masses in the centers 
of nearby gal axies and the properties of their hosting 
galax i es (e.g. Ferrarese & Merritt 2000; Graham et alj 
20011 iMarconi & Hunt ,2003; Gebhardt et al.. .20o6b 



have moved Active Galactic Nuclei (AGN) into the 
center of attention as key players in galaxy evolution. 

However, while it is clearly established that AGN 
must play a major role in galaxy formation to be able to 
produce the relation between supermassive black holes 
and their host galaxies at redshift zero, it is still a mys- 
tery when the correlation came into place and through 



which mechanism. To understand this, one has to un- 
derstand both the evolution of AGN and their connec- 
tion with their hosting galaxies. 

While our knowledge of low-redshift AGN and 
their correlation with the host is relatively broad, much 
less is known about high-redshift AGN hosts. Stud- 
ies of AGN hosting galaxies seldom reach beyond red- 
shifts of two. Resolving the host galaxy becomes more 
and more challenging as the faint hosts are extremely 
hard to disentangle fr om the nuclea r emission (e.g 



Hutchings 2003; Jahnke et alj l2004t ISchramm et al, 



2008l;IVillforth et alj|2008h . This is due to both the de 
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crease in apparent host galaxy size and surface bright- 
ness dimming. Low-luminosity AGN have a much 
more favorable core-to-galaxy ratio, making it possi- 
ble to study the properties of their hosting galaxies up 
to very high redshifts. 

Another open question is the general evolution of 
AGN over redshift. It is known that the luminosity 
function of high-redshift AGN significantly differs 
from the one at low redshifts. High-redshift AGN 
are on average much brighter than their low-redshift 
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counterparts dPunlop & Peacock However, 
given the fact that almost nothing is known about low- 
luminosity AGN at high-redshift, this is just the tip 
of the iceberg and we know little about the shape and 
normalization of the luminosity function on the faint 
end. 

New studies also imply that high-redshift AGN 
might be intrinsically different from their low-redshift 
counterparts. Jiang et al. (2010 ) found an interesting 
sample of dust-free AGN at high-re dshift that have 
no low-redshift counterparts. And IShemmer et al 



(b009) found a sample of weak-lined high-redshift 
AGN that do not seem to have low-redshifts counter- 
parts. On the other hand, the interesting AGN subclass 
of BL Lacs, highly variable objects thought to be the 
beamed counterparts of F anarof-Riley I radio galaxies 
(lUrrv & Padovani Il995h so far have only been de- 
tected at low-redshift. The missing high-redshift BL 
Lacs are st ill a mystery given their extreme intrinsic 
brightness ( Stockd 2001 '). Latest theoretical models 
for the evolution of AGN imply that the missing BL 
Lac problem might be due to the fact that certain types 
of AGN preferably appear at certain redshifts, with 
Fanarof-Riley t ype I galaxies being extremely rare at 
high redshifts ( Garofalo et al. 2010l) . Studying low- 
luminosity AGN at high-redshift might therefore also 
help answering the question what determines the in- 
trinsic properties of AGN. 

It is also still a puzzle how supermassive black holes 
formed. The fact that AGN at redshifts greater than six 
are fou nd to host superma ssive black holes as large as 
IO^Mq (' jiang et aL 2010 *) is rather puzzling given the 
fact that the universe was less than a gigayear old at 
those redshifts. This poses the question of how and 
when those black holes were formed. Getting a more 
complete view of black hole masses at high redshift 
might help understand how, when and in which objects 
mass accumulation took place. 

Despite the great interest in this topic, our knowl- 
edge about high-redshift AGN is still limited and al- 
most nothing is known about faint AGN at high red- 
shifts. To learn more about this topic, we need new 
complete samples of high-redshift AGN reaching to 
much lower luminosities. Deep multi-band surveys of- 
fer the possibility to select such samples. 

AGN can be identified in several ways, the most 
obvious being through spectroscopy. AGN can ei- 
ther be identified through extremely broad lines 
or by determining lin e ratios (Baldwin et al. 19811 
Veilleux & Osterbrockl Il987l) . The former has the 



problem that it only selects the small subset of broad- 
lined Type I again and misses Type II AGN as well 
as weak-lined objects. Both have the problem that 
spectroscopy is extremely costly especially for faint 
sources. Therefore, while spectroscopic surveys can 
play an important role in understanding high redshift 
AGN, carefully selected candidate samples are neces- 
sary to make this method efficient. 

A commonly used method for AGN selection is 
through their optic a l colo rs. This method was first 
used by iMarkarian (ll967^ who selected objects with 
excess UV flux and created the first catalog of nearby 
AGN. Possibly the most famous color selected AGN 
catalog is the Palomar Green Bright Quasar Catalog 
(Schmidt & G reen 1983, known under the name BQS 
or PG) in which objects with U - B < -0.44 mag were 
selected. However, UV excess might also be caused 
by star formation. Therefore more sophisticated meth- 
ods are n eeded to avoid contam ination by star forming 



galaxies. IWarren et al.l (119911) designed a multi-color 



selection method, in which only point-like sources are 
chosen and stars are rejected by their specific location 
in a multi-dimensional color space. Similar methods 
have been used for quasa r selection in the Slo an Dig- 
ital Sky Survey (SDSS) dRichards et al.ll2002l) . Such 
methods are observationally cheap and can easily be 
applied to big multi-waveband surveys. However, they 
rely on considerable deviation from a normal stellar- 
dominated spectral energy distribution (SED). There- 
fore such methods are not suitable to detect the in- 
teresting sample of low-luminosity AGN in which the 
galaxy emission dominates the overall SED. 

Other ways to select AGN are through excess radio. 
X-ray or mid-infrared emission. Radio selection has 
the downside that it only selects a small subsample of 
AGN. Only about 10% of all optical ly-selected AGN 
show considerable radio emission dSmith & Wright 
1980 !). Radio surveys are also generally shallow or 



may suff'er from confusion due to large apertures. An 
expection to this rule is the FIRST Survey which im- 
aged a large area of the sky using the NRAO Very 
Large Arry, the resolution of this survey is similar to 
ground-based optical surveys and the de pth of the sur- 
vey e xceeds most other radio surveys (^ Becker et al, 



19951) . X-ray surveys have similar problems to most 



radio surveys, the resolution is generally rather poor, 
causing problems with confusion. Due to the limited 
size of space-based X-ray telescopes it is also hard to 
detect fai nt high redshift sources (Chandra has a 1.2 
m mirror JWeisskopf et al.ll200(3l) . the XMM-Newton 
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effective collecting area corresponds to a mirror size 
of only about 80cm ('ja nsen etal ] l200lh l. 



Mid-infrared emission is thought to originate from 
warm dust in the obscuring dust torus surrounding the 
AGN (Sander s et al.lll989b . However, also star form- 
ing galaxies such as ultra luminous infrared galax- 
ies (ULIRGs) sho w excess mid-infrared emission 
(jGenzel et al.lll998l) . As star forming galaxies are es- 
pecially common at redshifts around 2 or 3 (see e.g. 



Madau et al. 1998: Reddv et al. 2006: Bouwens et al. 



2010), mid-infrared selected samples are polluted 
by star forming galaxies, especially at higher red- 
shifts. Additionally, at very high-redshifts, the Uni- 
verse might have been too young for considerable 
amounts of dust to be produced. Indeed, some high- 
reds hift AGN have be en shown to be virtually dust- 
free dJianget al.ll201ol) . 

Given that practically all AGN vary on all timescales 



from hours to decades (for a review, see lUlrich et al. 



1997h . variabi lity can be used as a selection criterion 



for AGN (e.g. iSaraiedini et al.l2003b . Due to light time 



travel arguments, any variability detected in galaxies 
on human-observable time-scales must originate from 
the nuclear region. Interestingly, it has also been found 



counterparts (e.g. 


Trevese et al.' 


'1994'; 'Cristiani et alj 


19961 di Clemente et al., 


,1996; ,de Vries et al. 
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Wold et al. 2007b. This makes variabihtv-selected 



samples especially sensitive to the interesting and oth- 
erwise difficult to detect sample of low-luminosity 
AGN. Additionally, it has been found th at AGN vary 



more strongly on shorter wavelength (e.g. d i Clemente et al. 



ri996) . This makes variability selection more sensitive 
at higher redshift in a given optical waveband. 

Deep multi-waveband surveys such as th e Great 



sampling a time-span of 5 years. They found evi- 
dence for nuclear variability in 16 galaxies down to 
a V-band magnitude of 27.5. Given that this study 
was performed on only two epochs of data, this was 
an extremely encouraging result, showing that vari- 
ability selection can succeed even for sparsely sam- 
pled data. The AGN found by lSarajedini et alJ (l2003h 
show redshifts up to 1.8. They were able to show that 
the variability selected AGN cover a wide range of 
colors and might not have been detected using color- 
color selection criteria. Based on their sample they de- 
rived an AGN luminosity function and found that low- 
luminosity AGN are possibly more abundant at high 
redshifts. 

Cohen et al. I (l2006b published a similar study based 



on the Hubble Ultra Deep Field (HUDF) /-band data, 
sampling a time-span of about 4 months. They found 
about 1 % of the sources to show significant variability 
with photometric redshifts as high as 5. 



Klesman & Saraiedinil (l2007h presented a study of 



the five epoch V-band data in the GOODS South field. 
They limited their study to infrared power-law and X- 
ray selected sources. As much as 26% of their AGN 
candidate sample showed variability. 



Trevese et al 



Observatories Origins Deep Survey (GOODS) (IGiavalisco et al 



20041) are very well suited for creating such a variabihty- 
selected AGN sample. Data for such surveys are typ- 
ically taken in several epochs, distributed over several 
months and are therefore suitable to detect variabitil- 
ity on timescales of months. Additionally, extremely 
deep imaging on a wide range of wavelengths all the 
way from X-ray to radio gives the possibility to study 
the broadband multiwa velength properties of variabil- 
ity selected AGN (e.g. iPaolillo et al.ll2004b as well as 
their parent population. 



Saraiedini et al.l (l2003b were the first to attempt as- 
sembling a variability selected sample from multi- 
epoch survey data in the GOODS Fields. They used 
the two-epoch Hubble Deep Field (HDF) V-band data. 



( 12008b applied similar method to 
the Southern inTermediate Redshift ESO Supernova 
Search (STRESS) data. This survey covers about five 
square degrees around the Chandra Deep Field South. 
In contrast to previous studies discussed here, the data 
were taken with a ground-based telescope. This results 
in significantly lower resolution, making it harder to 
detect the variability of a point source against the flux 
of its hosting galaxy. They applied their method to 
eight epochs of V band data taken around the Chandra 
Deep Field South (CDFS) and found 1 12 out of 5138 
.objects (about 2% of all sources) to be variable. 

Detecting variability from sparsely sampled data is 
extremely challenging. In general, there are two ap- 
proaches to variability detection. The first approach 
is to use well calibrated and robust statistical estima- 
tors such as the statistics. This approach allows for 
an effective control of false positive rates as the es- 
timators have known expected distributions for sam- 
ples of non-variable sources. The downside of such an 
approach is that one has to rely on the correctness of 
the error measurements. Given the fact that all kinds 
of unknown sources of error and systematics might 
influence aperture photometry, this approach is often 
avoided. The only statistical estimator for variability 
selection that does not show this drawback is the ANal- 
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ysis Of VAriance (ANOVA) dde Diegoll2()10l) . It can 
however only be applied to data with large numbers of 
data points that is oversampled in the time domain (see 
e.g. lVillforth et alfcoOOUde Diegoll20ld) . 

The second approach derives 'typical' scatter at a 
certain magnitude and then declares a certain number 
of objects that show 'significantly mo re scatter' vari- 
able (this approach has been taken by Saraiedini et alj 
2003; Cohen et al. 2006; Klesman & Sarajedini 20Q3 
Trevese et alj 20081) . While this approach is robust in 
a sense that the scatter in the 'real' data is taken into 
account, it is lacking a control of false positives and 
locally differing errors. 

In this paper, we will present a method that com- 
bines the advantages of both approaches. Our method- 
ology takes into account the observed scatter but uses 
well-calibrated statistical estimators. We are therefore 
in full control of expected false positive rates. 

We present a catalog of 139 variability-selected 
AGN from the GOODS North and South field in the 
Hubble Space Telescope (HST) F850LP band. This 
is the first variability-selected AGN catalog with a 
known expected contamination rate. Properties of the 
variability-selected AGN will be discussed in an up- 
coming paper 

Sample selection and data reduction are presented 
in Section 2. We present and discuss statistical meth- 
ods for variability selection in Section 3. We present 
the application of the discussed methods to our data set 
in Section 4, followed by Discussion in Section 5 and 
Summary and Conclusions in Section 6. 

2. Sample Selection and Data Reduction 

For our variability study we select all objects from 
the G OODS five epoch z-band catalog jGiavalisco et al. 
20041) with a signal-to-noise of 20 or greater. This 
results in a magnitude limit of about 25-26 in the 
Z/F850LP band. The z band was chosen as it has the 
deepest imaging of all the space-based optical observa- 
tions. Another factor in choosing the most appropriate 
waveband is the fact that AGN are on average more 



variab le at shorter wavelengths (e.g. ldi Clemente et al 



19961) . For this study, we decided to use the z-band as 
it is considerable deeper than the redder bands. Ad- 
ditionally, using blue wavebands might increase the 
influence of star forming regions, especially at higher 
redshifts and therefore diminish our detection power 
The signal-to-noise cut-off showed to be appropriate 
for this study as the percentage of variability-selected 



objects at magnitudes greater than 25 starts to be ex- 
tremely low. This indicates that the signal-to-noise at 
fainter magnitudes is too low to pick up variability. 

Data have been taken with the Advanced Cam- 
era for Surveys (ACS) Wide Field Channel (WFC) 
aboard Hubble Space Telescope (HST) in the F850LP 
band. Data reduction is performed using MultiDrizzle 
dKoekemoer et al.ll2002h . 

Photometry is performed using the NOAO-IRAli 
aperture photometry package daophot. Apertures with 
four different radii (0.12, 0.24, 0.36 and 0.72") cen- 
tered to the catalog positions are used. Different aper- 
ture sizes are used to determine the best aperture for 
this study. Several factors are important when se- 
lecting apertures, those are: signal-to-noise ratio, in- 
fluence of small changes in PSF and influences of 
galaxy light included when detecting variability of a 
point source against its hosting galaxies. The smallest 
aperture was chosen becuase at even lower radii, PSF 
changes become dominant, the largest was choosen be- 
cause initial tests showed the signal-to-noise ratio to 
plummet around the value choosen for the largest aper- 
ture. 

The flux in units of counts per second is derived 
for all five epochs and all four apertures used. This 
is done by summing over the aperture and subtracting 
the sky determined in an annulus with an inner radius 
of 3"and a width of 1 "using the NOAO-IRAF sky al- 
gorithm 'ofilter'. This algorithm uses the optimal fil- 
tering algorithm and a triangular weighting function 
employing the histogram of the sky pixels. It is found 
to yield to show the least bias in the presence of faint 
objects in the sky region (Henry C. Ferguson, private 
communication). 

Flux meas urement errors are de termined using 
MultiDrizzle dKoekemoer et al. I I2OO2 !) weight maps. 
Inverse variance weight maps are summed over the 
aperture. These measurement errors include errors in- 
troduced by the data reduction process. To check for 
the correctness of the weight maps, we perform the 
following test. Blank sky positions are selected and 
the flux in the blank positions are measured over five 
epochs. The errors for the sky positions are then de- 
termined from the weight maps. No variability is de- 
tected using this process. We therefore conclude that 
the weight maps reflect the errors from the data reduc- 
tion process correctly. If the weight maps would be 
incorrect, they would either over- or underestimate the 



~ http ://iraf . noao.edu/ 
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flux measurement errors. As there should be no vari- 
ability in blank sky positions, this procedure is suitable 
to determine if the measurement errors derived from 
the weight maps are correct. 

Background subtraction errors are included by 
quadratically adding the standard error of the mean in 
the sky ring times the number of pixels in the aperture 
to the flux error from the weight maps. To account for 
shot noise, exposure time corrections are determined. 
To do this, we examine the distributions of the flux 
maesurement errors from the weight maps derived for 
the objects. As the error determined from the weight 
maps only includes errors introduced by the data re- 
duction process, the flux error distribution shows three 
separate peaks, corresponding to areas with four, two 
or one pointings per epoch (i.e. longer exposure times 
correspond to lower errors from the data reduction 
process). According to the measurement errors deter- 
mined by the weight maps, the corresponding expo- 
sure time is assigned for each object in each epoch. 
Using this exposure time, we calculate the full stan- 
dard deviation and variance from shot noise. Standard 
deviation and variance are scaled to units of counts 
per second. The standard deviation is calculated as 
follows: 



are compar ed to methods previously used in simi- 



yjicr^ht X texp)^ + flux X t,_,p 

O-full = (1) 

Where cr,j,/„ is the measurement error from Mult- 
Drizzle weight map, flux is the flux in the given aper- 
ture and texp is the inferred net exposure time. 

The scaled variance is calculated as follows: 



varfuii 



(O" wht X texp)^ + flux X tgxp 



(2) 



exp 



Note that we calculate both the standard deviation 
and the variance scaled to units of counts per second. 
Due to the fact that the square root is not a linear opera- 
tor, var = (T^ is only valid in units of electrons. There- 
fore, when fluxes are scaled to counts per second, the 
corresponding standard deviation and variance have to 
be calculated separately as var = cr^ no longer holds. 

3. Statistical Methods for Variability Detection 

In this section, we present and discuss differ- 
ent statistical methods suitable for variability de- 
tection in sparsely sampled data. These methods 



20031: ICohen et al.ll2006t 



Trevese et al .1 120081) . Sta 



lar studies ("Saraiedini et al 
.Klesma n & Sarajedini 2007; 
tistical calculations are performed using the Python 
numerical and scientific packages NumP^andSciP>0. 
NumPy random number generators are used for all 
simulation^. 

For sparsely sampled data, only few statistical 
methods are suitable for variability detection. In this 
paper, we test and discuss the;^'^, F and C statistics. 

The statistics is widely used in general model 
fitting. In variability detection, it is used under the 
null-hypothesis that a flat line fits the data, which cor- 
responds to a null hypothesis that the object is not vari- 
able. It is calculated as follows: 



X 



(3) 



Where y is the flux, y is the mean over all y, and 
(Tj is the measurement error associated with a given 
flux measurement y,. This method has the advantage 
that it associates each flux measurement with its error 
estimate (as we will see soon, this is not the case for 
the other two methods). It is therefore well suited if 
one expects the errors to show significant deviations 
among epochs. 

The F statistics compares the expected to the ob- 
served variance. It is calculated as follows: 



^'^^observed 
VCirgxpecred 



(4) 



Where varghserved is the variance in the flux mea- 
surements and var expected is the mean of the flux error 
estimates given as variances. Note that, opposed to the 

statistics, the measurements errors are not associ- 
ated with individual flux measurements. Instead, both 
the flux measurements and error estimates are consid- 
ered as samples. The F statistics should therefore be 
restricted to cases in which similar errors for each mea- 
surement are expected or observed. 

The F statistics is also used in the ANalysis Of 
VAriance (ANOVA). This statistical method is a pow- 
erful tool for variability detection in cases in which 
the time domain is oversampled and many measure- 
ments are available. ANOVA does not rely on er- 



' http ://numpy. scipy.org/ 
http ://w w w. scipy.org/SciPy 
^Numpy random seeds from /dev/urandom 
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ror measurement but derives the expected variance 
from sub-samples of the data. It is therefore used 
in micro- yariabiHty studies (Ide Diego et alJ Il998 : 
Ivillforth et al.ir2009l) and is shown to have a very high 
power (where power is defined as o ne min us the rate 
of Type II errors or false negatives) (Ide Di ego 2010). 
ANOVA can however not be used in this study due to 
the small number of data points available. 

The C statistics compares the expected to the ob- 
served standard deviation. It is calculated as follows: 



C 



^observed 
^expected 



(5) 



In our case, (Tobserved is the observed standard devi- 
ation in the fluxes and CTe.xpected is the mean of the flux 
error estimates given as standard deviations. Note that 
this is mathematically not identical to the F statistics 
where the variances are compared. From a mathemat- 
ical standpoint it should be noted that opposed to the 
variance, the standard deviation is not a linear opera- 
tor. For example, when adding errors, variances can be 
added directly while standard deviations cannot. Just 
as the F statistics, this method does not relate each 
measurement to its error and should therefore only be 
used when similar errors are expected for each epoch. 

Previous studies of similar data sets have used other 



methods. ISarajedini et alJ (120031) analyzed two epochs 
of y-band data of the Hubble Deep Field (HDF). They 
calculated the difference in magnitudes between the 
two epochs and declared all objects with a deviation 
three times bigger than the standard deviation of all 
differences at a given magnitude as variable. They 
found the distribution of the differences to be close 
to Gaussian and therefore used 2-tailed p-values for 
a Gaussian distribution to estimate the number of false 
positives. This method seems very robust, it also pro- 
vides false positive estimates. 

Cohen et al.1 (l2006l) analyzed four epochs of data 



from the Hubble Ultra Deep Field (HUDF). They com- 
pared each of the measurements in the four epochs 
to the measurements in all other epochs, yielding six 
measurement pairs per object. They then empirically 
determined the error distribution for each epoch pair 
and declared all objects with a difference bigger than 
3cr variable for the given epoch pair Variability is then 
derived by analyzing the magnitude differences in each 
This i s com parable to the method used 
20031) . 



of the pairs 



by Saraie dini et al.l (120031) . However, they only had 
two epochs of data and therefore were limited in their 
choice of appropriate statistical methods. 



The method used by ICohen et al. (l2006h has the 
downside that variability is estimated from six differ- 
ent values and the ultimate selection criterion used 
is rather complex. Determining false positive rates 
would be possible using simulations, but is not straight 
forward due to the complex selection criterium used. 

A flavor of the C statistics was used by both 



Trevese et alJ (120081) and lKIesman & Sarajedinil (l2007h 



. In this method the standard deviation of the mea- 
sured magnitudes of an individual object is compared 
to the mean standard deviation of the entire sam- 
ple under the assumption that most objects are not 
variable. The two studies differ slightly in the way 
the y calculate the in ean error for each object. While 
Trevese et alj (2008|) average over the standard devia- 



tions, |IG£sm^n&_Saraiedini| ( l2007n parameterize vari 



ability using a quantity derived from the average of the 
variances in the magnitudes measured for each epoch. 
While it is true that var = cr^, it should be kept in mind 
that a square root is not a linear operator Therefore, 
var = cr^ is only valid in the space of the measure- 
ment. This equation is however no longer valid when 
the measurement is transposed into another space like 
counts per second (which represents a multiplicative 
transformation from the original space of the mea- 
surements) or magnitude (which is a highly non-linear 
transformation)!^. 

In those previously mentioned studies, sources that 
show a scatter of more than e.g. three times the normal 
scatter are then labeled variables. These methods are 
similar to C statistics under the assumption that ah ob- 
jects at a given magnitude have the same errors (i.e. the 
errors are fully shot-noise dominated and the net expo- 
sure times are identical for all objects and epochs). 

This method has several problems and caveats. Us- 
ing this method, it would require extensive simulations 
to associate the chosen selection criterion (in the exam- 
ple, 3cr) with a significance value. Such simulations 
have not been performed by the authors. Therefore, 
no expected false positive rates are avialable for those 
studies. This caveat can be important for large samples 
if they have a high number of false positives, in which 
case it would adversely impact the statistical signifi- 
cance of the number of detections. The statement that 
there are 100 variables technically only means that the 



*The reader can easily verify this by thinking about the variance and 
standard deviation of a flux measurements with 100 counts in 10 
seconds that needs to be transferred into counts/second space. It 
will become clear that var = o~ cannot be true in the transformed 
space. 
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100 most variable objects are listed. The fact that the 
variability limit is derived from the data also intro- 
duces the problem that more variable sources will re- 
sult in a higher limit, i.e. the detection limit depends 
on the number of variable sources. This will for ex- 
ample result in an overdetection of variability in faint 
sources, caused by the fact that this method enforces 
similar percentage detection rates at every magnitude. 

Additionally, using this method gives little control 
over the correctness of the error measurements. While 
the method takes into account the real scatter of the 
data, it cannot determine if the measured errors are cor- 
rect. Therefore, point spread function (PSF) changes 
or location dependent errors will remain unnoticed us- 
ing this method. Apprehensions concerning erroneous 
error estimates seem very realistic given the complex- 
ity of the errors that can be introduced in the data re- 
duction process and the known complexity of space- 
based PSFs. 

On a side-note, using magnitudes for variability de- 
tection seems to be very common. However, errors in 
magnitude space will be both non-Gaussian and asym- 
metric, while being close to Gaussian in flux space. 
When using common statistical estimators, this can 
cause problems as most statistical methods explicitly 
or in-explicitly assume Gaussian errors. 

3.1. Calibrating statistical methods 

To compare the different methods, the p-values for 
the different estimators need to be calibrated, p-values 
are values assigned to a given value of an estimator 
(in our case x^^P ™d C). The p-value for a given 
value of the estimator is the probability that a value 
as 'extreme' occurs from random data. Two different 
p-values exist, and they should not be confused. 

Two-tailed p-values give the probability that a value 
derived from random data is further away from the 
center of the distribution than the given value or val- 
ues. Two-tailed p-values are symmetric for symmetric 
distributions. For assymetric distributions they are de- 
fined such that each tail contains the same probability. 

One-tailed p-values give the probability that a value 
derived from random data is bigger than the given 
value. This is the p-value to be used for variability 
studies. It can be derived for all distributions, includ- 
ing asymmetric ones. 

Wheter one- or two-tailed p-values should be used 
depends on the variability statistics used. For distru- 
butions in which both tails represent extreme variabil- 



ity (for example, if deviations from catalogue magni- 
tudes are given), two-tailed p-values should be used. 
For distributions in which one tail represents extreme 
variability while the other tail represents low variabil- 
ity, the one-tailed p-values should be used (for example 
when using the;\f^ statistics where low values indicate 
low variability and high values represnt high variabil- 
ity). 

The p-value associated with the value at which ob- 
jects are considered variable is equal to the signifi- 
cance. The significance determines the number of false 
positives or Type I errors. It will also influence the 
power of the test. The power describes the percentage 
of real variables detected variable, and therefore is re- 
lated to the number of false negatives (Type II errors). 
Choosing a very strict significance limit will result in 
a low power On the other hand, choosing a very loose 
significance limit will result in a higher power. There- 
fore, when comparing different methods one has to be 
sure to compare them at the same significance. 

Some authors use 3cr or 5cr to describe the signifi- 
cance levels used in their studies, defining objects that 
deviate from the commonly observed scatter by more 
than 3/5o- as variable. This is consistent in itself but 
often wrongly related to p-values. For example, relat- 
ing 3cr to a significance level of 99.7% is wrong in the 
concept of variability detection as this is the two-tailed 
p-value for the Gaussian distribution. As discussed 
above, using two-tailed p-values is correct only for cer- 
tain statistics. And this still leaves aside the fact that 
the distribution is likely not Gaussian and therefore the 
association between the cr values and the p-values for 
Gaussians does not hold. 

The p-values for the three statistics used in this 
study at different significance levels are shown in Table 

m 

3.2. Testing the power of different statistical tests 

Using the p-values derived in the previous section, 
the power of the different statistical tests can be de- 



Table 1: One-tailed p-values for estimators F and 
C. P-values are given for light curves with 5 data 
points at different significance values. 



Estimator p = 95% p = 99% p = 99.9% p = 99.99% 





9.48 


13.26 


18.46 


23.41 


F 


1.89 


2.65 


3.69 


4.70 


C 


1.38 


1.63 


1.92 


2.17 
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termined. Mock variable light-curves are created, ran- 
domized and the statistical estimators are derived. De- 
tection rates are then calculated. Three different types 
of mock variable data sets are studied: 

• Noise: the flux for each data points is drawn 
from a normal distribution with a range of stan- 
dard deviations cr,,^,^ 

• Slope: the underlying theoretical light-curve is 
a simple linear trend with a range of slopes, each 
data point is randomized with a standard devia- 
tion of (To 

• Burst: the underlying light-curve is completely 
uniform with only a single outlier with a range 
of differences between base and peak value, 
each data point is randomized with a standard 
deviation of cTq 

A measurement error of cTq is assigned to all flux 
values. A range of 'variability strengths' Y is stud- 
ied for all mock lightcurve classes. The variability 
strengths V are defined as follows for the different 
mock light-curve classes: 



Noise: V : 



• Slope: y 



^var 
(To 

tnax(lightcurve) — min(lightcurve) 

(To 



• Burst: V = 



peakvalue — basevalue 



The power of the different statistical tests (i.e. the 
detection rate) at a given significance is then derived 
for the different types of light curves for a range of 
variability strengths. The detection rates for three 
methods show only minimal deviations that are con- 
sistent with error expected from the fact that we only 
use two decimals and therefore our p-values are not 
completely accurate. The results are shown in Table|2] 

3.3. Testing robustness of tests 

As we have shown, the three different statistical 
tests have equal power for 'perfect' data. However, 
real data is hardly perfect. Therefore, we will deter- 
mine how robust the statistical tests are in presence of 
deviations from the perfect simulated data. 

In the previously presented test, the flux measure- 
ment errors (meaning the standard deviation or vari- 
ance) are measured accurately. However, this is not 



expected for real data. Errors are measured in a similar 
way to fluxes and this process introduces measurement 
errors also into the error measurement. Therefore, it is 
of interest to see how different tests perform for 'er- 
roneous' error measurements. For our mock data that 
means that the flux measurement errors will also be 
drawn from a normal distribution with a certain width. 
The width of this distribution will give the 'defective- 
ness' of the error measurements. 

When rerunning the tests using 'erroneous' errors 
we find that the different tests differ in their detection 
power The statistics now shows the highest de- 
tection power, the C statistics shows the second high- 
est power and the F statistics shows the lowest power. 
Results for the detection power with 'erroneous' errors 
are shown in Table |4] located in the Appendix. 

We therefore conclude that the statistics should 
be used in cases in which the errors show strong devi- 
ations between the different epochs. The C statistics is 
to be preferred over the F statistics due to its greater 
power under the influence of erroneous error measure- 
ments. We will therefore only use the C and;^f^ statis- 
tics for our study. 

3.4. Estimating the influences of sparse sampling 

Additionally, it is of interest to understand the de- 
tection power for AGN light curves. To determine the 
detection power, we create mock AGN light-curves 
using the method introduced by iTimmer & Koenia 
(Il995b . This method randomizes both the phase 
and the amplitude of the Fourier transform. In other 
methods, only the phase is randomized and therefore 
only a subset o f all possible lightcurves is simulated 
dTimmer & Koenig, 199 5). 

10000 mock AGN light curves are created. We 
simulate light curves ten times longer than the sam- 
pling time-scale to in clude the red noise leak (see e.g. 
Vaughan et al.ll2003l for a discussion of the red noise 
leak), five data points are drawn from the mock light 
curves with sampling similar to the GOODS five epoch 
time sampling. The data are then randomized with a 
range of errors and data are analyzed using the three 
statistical estimator discussed. 

As a measure of the variability strength, we give 
the ratio between the assigned flux measurement errors 
and the median spread in the analyzed light curves in 
percent. With changes on a timescale of years of typ- 
ically around 1 mag in AGN light curves, a common 
error of 0.1 mag would result in a signalstrength of 
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Table 3 

Detection power (in per cent) for mock AGN light curves for different significance levels and different errors in 

per cent as defined in the text. 



eiTor in % 


Power (95%) 


Power (99%) 


Power (99.9%)) 


0.11 


100.00 


100.00 


100.00 


0.23 


100.00 


99.97 


99.97 


0.34 


99.97 


99.97 


99.95 


0.45 


99.95 


99.95 


99.93 


0.57 


99.93 


99.92 


99.83 


0.68 


99.92 


99.82 


99.52 


0.79 


99.82 


99.53 


99.29 


0.91 


99.56 


99.32 


98.87 


1.02 


99.44 


99.02 


98.25 


1.13 


99.22 


98.73 


97.40 


2.27 


91.99 


87.42 


81.46 


3.40 


78.12 


69.34 


60.51 


4.54 


63.25 


53.08 


42.69 


5.67 


49.47 


38.72 


28.92 


6.81 


37.74 


27.84 


18.82 


7.94 


28.81 


19.34 


11.81 


9.08 


21.23 


13.33 


7.32 


10.21 


15.58 


8.86 


4.11 


11.35 


10.92 


5.49 


2.13 


12.48 


7.80 


3.32 


1.21 


13.62 


5.24 


1.93 


0.53 


14.75 


3.36 


1.18 


0.28 


15.89 


2.09 


0.61 


0.13 


17.02 


1.43 


0.35 


0.03 


18.16 


0.88 


0.20 


0.00 


19.29 


0.43 


0.07 


0.00 


20.43 


0.33 


0.02 


0.00 


21.56 


0.19 


0.00 


0.00 



about 3, resulting in a detection probability of about 
78% at 95% significance. Detection probabilities for 
other typical errors and variability strengths can be de- 
rived from Table [3] 

4. Results 

4.1. Zero-point Calibration 

First, we check for possible zero-point offset be- 
tween epochs. We analyze offset in the zero-point for 
bright, point-like objects. Zero-point offsets are de- 
rived for aperture corrections (i.e. the correction ap- 
plied to determine the entire flux of the object instead 
of the flux in the aperture) and inter-epoch zero-point 
drifts. This is done by calculating the mean offset 
between the measured magnitude and catalogue mag- 
nitude for bright, point-like objects. The zero-point 
drifts for the large apertures are very mild, indicating 
changes below 0.01%. This implies that the photomet- 
ric calibration is extremely accurate, the differences 
are within the error of the estimator used. 

For the smallest aperture (0.12"radius) however, 
the fitting indicates inter-epoch changes on the level of 



~ 1%. We therefore correct the fluxes and errors with 
the derived zero-point offsets for the 0. 12"aperture and 
compare the number of detections with and without the 
correction applied. The number of variables increases 
when correcting for zero-point drifts indicating that ei- 
ther the error of the fitting routine is dominating or that 
changes in the PSF are location-dependent. Therefore, 
we test for possible problems with PSF instabilities in 
the next subsection. 

While we could simple use the biggest aperture to 
avoid such problems, this is not the optimal solutions. 
Big apertures can result in both low signal-to-noise ra- 
tios and lower detection rates for faint AGN due to 
high contribution from the hosting galaxy. Therefore, 
we will try to assess which aperture is optimal for this 
study. 

4.2. PSF Stability 

Using TinyTim ( lKristlll995h . we test possible influ- 



ences of defocusing over the field of view and changes 
in the PSF shape to get an estimate of the errors ex- 
pected from PSF changes. 

First, images are checked for possible changes in 
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the apparent position between the different epochs. 
Objects show decentering between different epochs 
only on sub-pixel scales, as this is comparable to the 
accuracy of the centering algorithm, no decentering is 
assumed between epochs. 

To test for possible errors due to changes in the PSF, 
PSFs for 64 ACS WFC positions with filter F850LP 
are created. Defocusing is introduced with values be- 
tween -5 and -1-5 fiin in steps of 1. These values de- 
scribe the movement of the secondary mirror with re- 
spect to the primary, where corresponds to the tele- 
scope being in focus. Those values are typical fo- 
cus changes due to 'breathing' of the HST spacecraft 
(Idi Nino et al. I l2008 h. Both a E galaxy template and a 
QSO template are used for the objects spectral shape, 
but no difference show in the resulting PSFs. Aper- 
ture photometry is then performed on each image with 
the same aperture sizes used for the data. From this 
measurement, typical errors due to both changes of the 
position of the object on the chip and the focus can be 
estimated. For the smallest aperture (0.12"), typical 
values for PSF changes expected are ~ 1%, dropping 
to ~ 0.25% for the 0.24"aperture, ~ 0.15% for the 
0.36"aperture and ~ 0.1% for the 0.72"aperture. 

Note that there is a known red halo in the ACS WFC 
for t he F850LP filter wh ich is not included in Tiny- 
Tim dSirianni et al .l2005h . This might slightly alter our 
results, making the possible changes due to defocus- 
ing and changes over the field of view smaller. False 
variability due to the changes in the red halo are only 
expected if spectral changes occur. In case the spec- 
tral shape changes, the object is intrinsically variable. 
Therefore, the red halo might cause true variability to 
be boosted, but no excess false variability is expected 
to be introduced due to this effect. 

Given the expected errors due to PSF changes, we 
will first asses if PSF instabilities induce false vari- 
ability for small apertures in the data. Figure [T] shows 
the detection rates for point-like and extended objects 
over the aperture radius used. As we can see the de- 
tection rate for point-like objects is generally signifi- 
cantly larger than for extended objects. This indicates 
that point-like objects are intrinsically more variable 
than extended objects. This is expected given the fact 
that most high-luminosity AGN would appear point- 
hke due to the fact that the AGN significantly out- 
shines the galaxy. 

However, the detection rate for stars is significantly 
higher at the two smallest apertures, indicating that 
PSF changes induce false variability. Figure |2] shows 



Table 2: Detection power (in per cent) for different 
mock light curves at a significance level of 99.9 per 
cent. V is the strength of the varaibility as defined in 
the text. 



V 


Slope 


Burst 


Noise 


0.5 


0.11 


0.17 


0.00 


1.0 


0.20 


0.35 


0.14 


1.5 


0.40 


0.97 


8.33 


2.0 


0.82 


2.46 


32.96 


2.5 


1.59 


5.74 


56.55 


3.0 


2.94 


12.14 


72.53 


3.5 


5.57 


22.62 


82.59 


4.0 


9.65 


37.15 


88.45 


4.5 


15.57 


53.80 


92.12 


5.0 


23.64 


69.91 


94.70 


5.5 


33.32 


82.97 


96.32 


6.0 


45.39 


91.74 


97.23 


6.5 


56.99 


96.60 


97.93 


7.0 


68.56 


98.81 


98.40 


7.5 


78.35 


99.62 


98.83 


8.0 


85.97 


99.92 


99.05 


8.5 


91.90 


99.98 


99.28 


9.0 


95.63 


99.99 


99.38 


9.5 


97.88 


100.00 


99.54 
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Fig. 1 . — Percentage of objects detected variable, sep- 
arated into point-like (red pentagons) and extended 
(blue circles) objects. Detections are performed using 
C-Statistics for four different aperture radii. 
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a comparison between the observed estimator distribu- 
tions of point-like and extended sources. The shape 
of the distribution differs strongly for the two small- 
est apertures. For the two biggest apertures however, 
the distributions agree well. This indicates that PSF 
changes will not affect those apertures. 

On the other hand, for the purpose of studying vari- 
able AGN, it should be noted that making the aper- 
ture bigger will start to drown the variability in the 
host galaxy light. This is due to the fact that bigger 
apertures will include more emission from the host- 
ing galaxy. Therefore, the smallest possible aperture 
should be used. As the two smallest apertures show 
signs of variability due to PSF errors, we use an aper- 
ture with a radius of 0.36 ". 

4.3. Re-calibration of statistical estimators 

The photometric errors of ACS have been studied 
excessively (Sirianni et al. 2005). However, given the 
extreme importance of exact error measurements in 
this study and the general complexity of photometric 
errors, we will check and if necessary correct the er- 
ror measurements derived. To check the quality of our 
measurement errors (cry^,,// as defined in the text), we 
compare the distributions of measured and C to the 
theoretical distributions. It is assumed that most ob- 
jects are not variable and therefore the observed distri- 
bution should follow the theoretical distribution. 

We find that the observed distributions for our mea- 
surements with an aperture radius of 0. 12 "agree rather 
well with the theoretical distributions (Fig. O, even 
given the PSF problems. Using bigger apertures, the 
distributions start to deviate from the theoretical ones, 
showing a flatter tail towards higher values than ex- 
pected (Fig. O. This indicates that either the shot 
noise or the background contribution are underesti- 
mated. This will result in excess false positives. How- 
ever, as we have discussed in the last paragraph, aper- 
ture radii of 0.36"should be used to avoid false posi- 
tives due to PSF changes. 

To check more closely for possible problems in the 
error determination, we compare the variability esti- 
mator distributions to the theoretically expected distri- 
butions for different object magnitudes. We find that 
the distributions shift to higher values at higher fluxes 
for all apertures (see Fig. |4]for an aperture radius of 
0.12"and left panel of Fig. |5]for an aperture radius 
of 0.36"). This could be related to an underestima- 
tion of shot noise errors or other factors. To correct 



for this problem, which could cause an over-detection 
of variability in bright sources, we decide to include a 
flux-dependent factor in the variability estimators. To 
derive this factor, we use the C statistics. When ap- 
plying a correction to the estimator, the error estimates 
are changed for all five epochs. This means that a cor- 
rection is applied to the sample of errors. This is more 
compatible with the C than the;^'^ statistics. Therefore, 
from now on the C statistics will be used. 

To derive the correction factor, we divide the sam- 
ple into seven flux bins. A histogram of the test statis- 
tics is then calculated for each bin and a gaussian is 
fitted to determine the peak position of the distribu- 
tion. This estimate for the distribution peak is then 
plotted against both the mean and the median of fluxes 
in each bin (see Fig. |6]l. A clear trend for the peak to 
move to higher values at higher fluxes is visible. As 
the increase is presumably due to an underestimated 
flux error, a square root function with an y-axes off- 
set is fitted to the data. The correction is then applied 
to the data by normalizing the test statistics with the 
value derived. From the definition of the C-Statistics, 
we can see that this can also be interpreted as a correc- 
tion applied to the error estimates. 

The flux corrected C-distributions are shown in the 
right panel of Fig. |5] As we can see the distribu- 
tions now agree well with the theoretical distributions. 
There are small deviations, but those indicate that the 
correction produced a slight overshoot, resulting in an 
overcorrection for very low C values. This could pos- 
sibly cause an increase in the number of false nega- 
tives. However, as the deviations are only apparent in 
the left tail of the distribution, the possible effects of 
this overcorrection should be minimal. Using the flux- 
corrected 0.36"C value gives us 173 variable sources 
out of 11931 total sources at a significance level of 
99.9%. 

4.4. Clustering of variable sources: testing for lo- 
cally underestimated errors 

Next, we will determine if location dependent er- 
rors or strongly location dependent changes in the PSF 
have been underestimated. Extreme clustering of vari- 
able sources might indicate such problems. These kind 
of problems could be very hard to pick up using the 
previously described test as they would potentially af- 
fect only a small percentage of sources. 

The spatial distribution of the variable objects in 
the GOODS North and South field is shown in Fig- 
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ures I2] and [8] At first visual inspection, a clear clus- 
tering at a few distinct locations is visible, especially 
in the GOODS North field. We inspect these areas by 
eye and find that those objects are located very close 
to bright stars that show clear diffraction spikes. The 
light curves show heavy outliers in some epochs. In- 
deed, visual inspection shows that diffraction spikes 
are present close to the centers of these objects. We re- 
ject the 18 objects affected by this problem, dropping 
our number of variables from 173 to 155. 

However, there might still be more subtle cluster- 
ing caused by underestimated location-dependent er- 
rors that is not picked up by eyeball inspection. There- 
fore, we calculate the distance to the nearest variable 
neighbor for all variables in the North and South field. 
We then create a mock variable catalog by randomly 
selecting the same number of sources for each field. 
Only objects with detections in five epochs are used 
for the selection. The distance to the nearest neighbor 
is also calculated for the mock catalogs. 

The distributions of nearest neighbor distances for 
the variable and mock variable catalogs are then com- 
pared using a two-sample Kolmogorov-Smirnov test. 
In case the variable sources would be extensively clus- 
tered, the p-values should be small, indicating that the 
samples are not drawn from the same parent popula- 
tion. 

We create 100000 mock catalogs for each of the 
fields. Kolmogorov-Smirnov tests are then performed 
for the variable sample against each of the mock sam- 
ples and against a master distribution created by merg- 
ing all mock distributions. The p-values for the North 
field range from 0.02 to > 99.9 per cent with a mean 
of 36 per cent and a median of 30 per cent. Comparing 
the distribution for the variables to the master distri- 
bution gives a p-value of 28 per cent. The p-values 
for the South field range from 10^"* to > 99.9 per cent 
with a mean of 32 per cent and a median of 28 per 
cent. Comparison with the master distribution gives 
a p-value of 20 per cent. This shows that no excess 
clustering is found after rejecting objects affected by 
diffraction spikes. This finding indicates that location 
dependent errors are well accounted for 

Auto- or 2-point-correlation functions might be 
more appropriate to study possible clustering, but as 
our test with nearest neighbor distances does not show 
any abnormalities, this easy test should be enough to 
check for problems with location-dependent errors or 
PSF changes. Additionally, weak 'real' clustering in 
our variability sample might be present. We will there- 



fore not further explore this topic. 

4.5. Rejection of saturated objects, stars and su- 
pernovae 

At this point, our catalog is only a catalog of vari- 
able objects, we are however interested in variable 
AGN. To produce a catalog of variable AGN, we need 
to reject saturated objects, stars and supernovae. 

Fig. |9]shows the location of the variable sources in 
a magnitude-stellarity plot. As a measure of stellarity, 
we use the Sextractor pa rameter ClassStar prov ided in 



the GOODS catalogue ([Dickinson et al.ll2003b . Ob 



jects with high stellarity can have two potential prob- 
lems. They could be saturated or the non-linear range 
of the chip could be reached, in both cases, false vari- 
ability is expected. 

To flag possible saturated objects or those in the 
non-linear regime, we determine at which magnitude 
the objects start to enter the non-linear regime. There- 
fore, the object flux is compared to the objects peak 
flux for all point-like sources and the turn-over point is 
marked. This happens at 18th magnitude. All objects 
showing a peak flux higher than the turnover point or 
are brighter than 18th magnitude and point-like are 
therefore rejected from our sample. None of our de- 
tected variable objects is in this regime. All bright stars 
were rejected before due to extremely high error bars, 
caused by the saturation or diffraction spikes. 

As for excluding stars, we have to take into ac- 
count that both stars and AGN can appear point-like. 
High-luminosity AGN can be more than 100 times 
brighter than their hosting galaxy. Therefore, they ap- 
pear point-like even if the host galaxy could theoreti- 
cally be resolved. Therefore, we will not reject point- 
like objects per se. 

However, we correlate our data set with other cata- 
logs to exclude stars. 15 objects in our variable cata- 
log turned out to be red stars and are therefore rejected 
from the variable AGN catalog. They are given in the 
final table but are flagged as stars. More of our vari- 
ability selected objects might be stars, but we are con- 
fident that we were able to flag the majority of stars. 

Supernovae are the only other objects that show 
variability on the sampled timescale. Therefore, we 
correlated our data with supernovae identified by 



Riess et al.l (120041) . This study identified supernovae 
from the five ep och GOODS data . We found that one 
supernova from lRiess et al.l (120041) is identified in this 
study as variable. This supernova (2003XX) went off 
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in the very center of an elliptical galaxy. The light 
curve of this object indeed shows a singly outlier This 
object is therefore rejected from the final variable AGN 
catalog. It is however listed in the variable table and 
flagged. 

4.6. Variable AGN catalog 

Now that we have determined a reliable variability 
estimator and rejected objects influenced by diffraction 
spikes, saturated objects, stars and supernovae, we can 
go ahead and asses the properties of the final variable 
AGN sample. 

Figure[TO]shows the flux-corrected C for a 0.36"aperture 
versus the object magnitude for the final sample. Vari- 
ability is detected down to magnitudes as faint as 25.5, 
with the faintest object being clearly detected with a 
significance of 99.99%. All of the variable sources 
brighter than 18th magnitude are galaxies. At fainter 
magnitudes, both point-like and extended objects show 
in our variabiUty sample, the most variable objects 
tend to be point-like. 

To assess if the signal to noise limit used in our data 
is too strict, the detection rate for the faintest sources 
is derived. For magnitudes > 25, the detection rate is 
close to the expected false positive rate. This indicates 
that a lower signal to noise limit would only result in 
more false positives and not more 'real' detections. 

With 11931 objects with five epoch detections in 
our catalog, false positive detections might pose a se- 
rious problem. To get an estimate for the expected 
false positives contamination, we conservatively as- 
sume that all objects are non-variable and calculate the 
number of false positives at a given significance. 

For a significance of 95% we expect as many as 596 
false positives. Even at a significance at 99%, the num- 
ber of expected false positives is still very high (119). 
Only for a significance of 99.9% does the number of 
expected false positives drop to 12. At a significance 
level of 99.99%, we expect only a single false positive, 
making a catalog with such a strict selection criteria 
'clean' from false positives. 

To see how this affects real data, we show the num- 
ber of raw and false positive corrected detections in 
a number of magnitude bins (Fig. fTTI) . False positive 
corrections are applied by assuming that all objects are 
non-variable and subtracting the expected number of 
false positives from the number of detections. We see 
that false positives pose a very serious problem at all 
significance levels lower than 99.9%. Thus, this very 



strict limit should be used when determining the vari- 
able object catalog. 

However, lower significance levels can be used 
when trying to estimate the number of variable objects 
in the field for statistical arguments. If one is not in- 
terested in knowing the individual objects that show 
variability but just the number of variable object in a 
given subsample, slack significance levels can be use- 
ful as they reduce the number of false negatives. This 
becomes evident when looking at Fig. [TT] there is a 
big number of false negatives at a 99.9% significance 
level. 

To estimate the true number of variables, we use a 
relaxed significance level of 95%, this gives us 1072 
'variable' objects, given the expected number of false 
positive (596) at this significance, we are left with 476 
'true' variable objects, from those, we subtract the 18 
objects found to show false variability due to spikes 
from bright objects. According to our findings, about 
3.8% of all objects are variable. This however, does 
not account for contamination due to stars. In our final 
99.9% significance catalog with 155 entries, 15 objects 
are stars. More objects might turn out to be stars, but 
could not be identified as such. Therefore, at least 10% 
of all variable objects identified in this study turn out to 
be stars. Assuming that 10% of all variable objects are 
stars, we derive that about 3.4% of all objects in our 
sample are variable AGN. (Note that the precentage 
of stars might depend on the magnitude and therefore 
our assumption of a fixed rate of stars might not be 
correct.) 

For the catalog, a significance level of 99.9% will 
be used. The final sample of variables therefore con- 
tains 155 objects, with catalog z magnitudes between 
16.45 and 25.5. As mentioned above, the conserva- 
tively estimated number of false positives is 12, reso- 
lution in an expected catalog contamination of about 
7.7%. Out of these 155 variable objects, 15 are identi- 
fied as stars and one object is a supernova, leaving 139 
variable AGN. 

Additionally, a 'clean' sub-catalog with a selection 
criterion of 99.99% is provided. This catalog contains 
93 objects and has only one expected false positive en- 
try. This sub-catalog contains 10 stars and one super- 
novae, leaving 82 AGN. 

Both final catalogs are shown in Table |5] A flag in 
the table indicated if the objects belongs to the 'nor- 
mal' (99.9% significance) or 'clean' (99.99% signifi- 
cance) catalog. The entire catalog will be made avail- 
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able at Viziefl 

4.7. Variability strength 

Additionally, we derive the intrinsic variability for 
the variability selected sample. The observed variabil- 
ity can be written as: 



^observed 



2 + 0-2, 



(6) 



Where cr is the measurement error. This gives the per- 
centage variability V: 



^intrinsic 

(flux) 



observed 



-<cr)2 



{flux} 



(7) 



The percentage variability V for all objects in our 
variable sample is shown in Figure[T2] Naturally, there 
is a lower envelope to the variability strength that is ris- 
ing to lower luminosity objects, caused by the fact that 
the signal to noise ratio is worse for lower luminosities, 
making it impossible to detect very subtle variability. 
V for all variable objects is also shown in the catalog 
table (Tables. 

To asses the detection limit, we fit a lower envelope 
to the data. The data is well fit by the following third 
order polynomial: 



V = 0.036xmag^-l.504xmag^+l5.993xmag-l.644 

(8) 

This equation holds for magnitudes > 21, for lower 
magnitudes, a lower detection limit of about 1 % vari- 
ability strength is found. 

5. Discussion 

5.1. Properties of variability selected AGN candi- 
dates 

All variability selected obj ects are matched to 
the GOODS s pectroscopic data ( !Vanzella etaP 12008 : 
Popesso et alj 12009,) . Redshifts, flags and absolute 



magnitudes of those objects are included in Table|5] 28 
out of 155 objects can be matched to the spectroscopic 
data. Out of those 28 objects, 5 have been identified 
as stars. For 3 objects, no spectroscopic redshift could 
be derived. The remaining 20 objects have redshifts 



between 0.045 and 3.7. Amongst those are five broad- 
hned AGN, with redshifts of 0.74,0.84,1.23,1.61 and 
2.80. With absolute rest-frame z band magnitudes 
ranging between -22.57 and -24.31, all spectroscop- 
ically identified AGN in our sample are rather faint. 
One spectroscopically identified object has an inferred 
absolute magnitude of -12.8, indicating that the red- 
shift identification might be faulty. However the red- 
shift estimation is based on a single line and might 
therefore not be correct. Practically all sources show 
extremely strong line emission, indicating either high 
star formation rates or AGN activity. 

All variability selected objects are visually in- 
spected. 41 objects are unresolved, out of those 15 are 
identified as stars, leaving 26 unresolved AGN candi- 
dates. 18 objects have a dominant core and show faint 
extended emission. Nine objects show clear signs of 
interaction, either in the merger stage or showing tidal 
tails. 16 objects are elliptical galaxies and four show 
clear disk structure. 41 are resolved, but no morphol- 
ogy can be determined. Those are mostly faint objects 
of small size. On e object (J033241. 87-2 7465 1.1) has 
been identified by Straughn et al. ('2006') as a tadpole 
galaxies. Those interesting objects are believed to be 



galax ies in an early stage of merging dStraughn et al, 
2OO69. However, this particular object is a rather ex- 



^http://vizier.u-strasbg.fr/viz-bin/VizieR 



treme example of this object class as it shows much 
less substructure than most tadpole galaxies. It looks 
similar to a strong lensed galaxy. However, there are 
no clusters nearby. 

Finally, seven objects show complex structures. 
Those galaxies either have multiple cores or are clearly 
extended with no clear center Note that objects with 
multiple centers might be wrongly identified as vari- 
ables due to PSF changes. This can happen even in 
cases in which the aperture is big enough to avoid 
false variability due to PSF changes for single-center 
objects. However, one of those complex sources 
(J033228.30-274403.6) is detected in X-rays. This 
object is also a B-band drop-out, indicating a redshifts 
around 2-3. Given that hot cluster gas emission is not 
detectable at such high redshifts, this indicates that 
there is indeed an AGN in this object! Therefore, 
complex objects will be included in the final catalog. 

5.2. Comparison witli otlier variable AGN cata- 
logs in the GOODS field 

Now that we have presented our variability-selected 
AGN catalog and the estimated occurrence of variabil- 
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ity, we would like to compare our sample to other sam- 
ples derived in a similar way in the same field. 



Sarajedini et al.l (i2003) studied HDF data and found 



16 variability selected AGN, none of those is found to 
be variable in our sample. They performed their study 
in the V-band, which should pick up diff'erent types 



of var iability than our z-band study ddi Clemente et al 



19961) . Additionally, they are covering a time-span 



of five years while our observations span only about 
a year. This will result in higher sens itivity for the 



detection of more lumin ous AGN (e.g . Treves e et al 
[L994; Wold et al. 2007). ISarajedini et a l. (2003) used 



the HDF data, therefore the field is much smaller than 
ours but the data is deeper All in all, their data is very 
diff'erent from ours, therefore, the missing overlap is 
somewhat expected. 

Cohen et al.l (|2006') studied /-band HUDF data and 
fo und 45 'best variabl e candidates'. The /-band used 
by Cohen et al. ( 20061) is only slightly bluer than the 
z-band used in this study. Therefore, the wavelength 
dependence of AGN variability between those bands 
should be negligible. However, they covered a time- 
span of only about four months, opposed to about a 
year in this study. This makes their study less sensitive 
to variability on time scales of months and longer. The 
area covered by ICohen et al. I (l2006h is much smaller 
than ours, but their data is significantly deeper Out 
of the 45 variable objects identified by ICohen et"al 



two are found in our catal og as variable . How- 
ever, the objects found variable bv lCohen et al. I (l2006l) 
are mostly fainter than 25th magnitude and the two ob- 
jects that are found in both catalogs are amongst the 
brightest in the Cohen catalog. Given that their data is 
much deeper, it is expected that little overlap is found 
between their and our sample. 



Klesman & Saraiedinil ( 120071) studied a sample of 
112 X-ray and infrared selected AGN candidates 
they selected from the GOODS V band five epoch 
catalogue. 29 of those objects showed variability. 
They used the V-band data as AGN variability is 
found to be stronger a t bluer wavelengths (see e.g. 



di Clemente et al.1 \1996l) . On the other hand, the 



band data used i n this study is deeper than the V- 
band data used bv lKlesman & Sa rajedini (2007). The 
time coverage for the iKlesman & Sarajedini (.20071) 
is practically identical to the one in this study. Their 
sample was drawn just from the GOODS South field 
and was also restricted to AGN candidates that were 
preselected, thus also accounting for the difference in 
sample size. 



Out of the 29 variables from the lKlesman & Sarajedini 

(E007 ) sample, 13 are detected variable in our study, 
12 of those are in our clean catalog, two objects are 
too faint to be included in our sample and two objects 
have been rejected due to suspiciously high error bars. 
Two objects have been detected variable in our study 
but not bv lkiesman & Saraiedinil ( |2007|) . All other ob- 
jects are not detected variable in both studies. This is 
a rather promising overlap. 

The differences might be due to several rea- 
sons. The different waveband used might play a role 
since variability is stron ger at shorter wavelengths 



( di Clemente et al. 19961). Additionally, we used a 



larger aperture than i Klesman & Sarajedinil ( 120071) . 
which we chose in order to reduce PSF-related vari- 
ability to less than 0. 15 percent. Finally, the variability 
detection for our z-band selected sample becomes lim- 
ited beyond about z 25, although a direct comparison 
with the V-band sample would also need to include 
differences in variability amplitude and color terms 
between the two bands. 



Trevese et al 



(I2008h selected variable AGN from 
the two year Southern inTermediate Redshift ESO Su- 
pernova Search (STRESS) data. Data were taken about 
every thre e months ove r a tim espan of two years. 
Therefore, iTrevese et al.l (120081) not only cover dou- 
ble the time span covered in this study, they also have 
more data points available, making it more likely to 
detect variability. They used V band data taken with 
the ground-based ESO/MPl 2.2 m telescope at ESO, 
La Silla (Chile). With a seeing of about 1", their 
resolution is about 10 times poorer than that in all 



other studies discussed here (ISarajedini et al 



Klesman & Saraiedinil 120071 : ICohen et al 



20061 



2003 



this 



study) . The poorer resolution makes ITrevese et al. 
(l2008h less sensitive to low luminosity AGN with 
clearly extended host galaxies. The area covered in 
their study is about six times larger than one GOODS 
Field and therefore about three times larger than our 
entire field. 



From the 112 variable objects found by Trevese et al. 



(12008) only 17 are also in the GOODS footprint. Out 
of those, three have been detected variable in our study, 
six objects did not have detections in all epochs, five 
are not variable in our study, two objects were too faint 
for our signal to noise cutoff and one object would 
have been labeled saturated in our study. Given the big 
differences between their and our study, such a small 
overlap is expected. 
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6. Summary and Conclusions 

In this paper, we have presented and discussed dif- 
ferent statistical methods that can be used to detect 
variable sources in large samples. Results have been 
applied to the GOODS North and South five epoch z- 
band data. 

All three statistical methods tested {C,F,x^) show 
equal powers for mock data, however, in cases in 
which the error determinations are erroneous, 
shows the highest power, followed by C and F. We 
have presented a robust statistical method for detect- 
ing variable objects from sparsely sampled data. Our 
method makes it possible to control the number of 
false positives and to test and correct for possible prob- 
lems in the error determination. 

We present a final catalog of 155 variable objects, 
selected with a 99.9% significance value. Out of those 
155 objects, 15 are identified as stars and one objects 
is a supernova, leaving 139 variability selected AGN 
candidates. The AGN candidates have magnitudes be- 
tween 16.5 and 25.5 mag in the z band. The cata- 
log has a expected false positive contamination of 12 
objects, corresponding to a contamination of 7.7% in 
the entire variable catalog. Additionally, we present 
a 'clean catalog' with a significance of 99.99%. This 
catalog contains 93 objects and 82 AGN candidates. 
It is expected to contain only a single false positive, 
corresponding to a contamination rate of about 1%. 
Detection rates at lower significance levels indicate 
that in total about 3.5% of objects are variable AGN. 
This is higher than the variables rates claimed in previ- 
ous studies (Saraiedini et al. 2003: Cohen et al. 20061 
Elesman & Sarajedinill2007l) . 

20 variability selected AGN have spectroscopic red- 
shifts. The redshifts of those objects lie between 0.045 
and 3.7. Amongst the 20 objects with optical spectra 
are five broad-lined AGN, with redshifts between 0.7 
and 2.8. With absolute magnitudes ranging between - 
22.57 and -24.31, all spectroscopically identified AGN 
in our sample are rather faint. This shows that this 
method is indeed suitable for detecting low luminosity 
AGN up to very high redshifts. 

Compared to other variability selected catalogs 
published so far, our GOODS catalog is larg er than 



the one published by ISarajedini et al.l (120031) on the 



Hubble Deep Field which covers a much s maller area. 



but go es deeper. Similarly, compared to .Cohen et al 



(120061) who covered the Hubble Ultra Deep Field, our 
sample also covers a larger area with about 3 times 



more sources, but is less deep. Since we study both 
GOODS fields and don't limit our sample to a pres- 
elected AGN catalog, our sample is also larger than 
that of iKlesman & Saraiedinil (|2007|) . and our cho- 
sen aperture size reduces our susceptibility to false 
positives. Our a rea co vered is smaller than the one 
in iTrevese et al.i (I2OO8I) . and they also have a better 
time sampling. However, their ground-based data has 
poorer resolution, making it much less sensitive to the 
interesting sample of low-lu minosity AGN. Finally , 
our study, along with that of ^Saraie dini et al (2003), 
are the only ones that provide estimates for false posi- 
tive rates. 

We have therefore presented a new extensive cata- 
log of variabiUty selected AGN in the GOODS Fields. 
Our catalog contains interesting low-luminosity AGN 
that cannot be detected using commonly used selec- 
tion techniques. This makes our catalog an interesting 
probe for the poorly understood population of high- 
redshift low-luminosity AGN. 

Given the formidable multi-wavelength coverage of 
both GOODS fields and the availability of both pho- 
tometric and spectroscopic redshifts, this sample can 
give unique insights into the properties of AGN up to 
high redshifts and open the window to a sample of pre- 
viously undetected low-luminosity high redshift AGN. 

The multi-wavelength properties and parent popu- 
lation of our variability selected AGN sample will be 
presented in an upcoming paper. 
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Fig. 2. — Observed distribution for variability es- 
timator C for point-like (red dashed line) and ex- 
tended (green dash-dotted line) objects. Objects are 
divided into classes according to their stellarity mea- 
sure ClassStar at a value of 0.5. 
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Fig. 3. — Comparison between observed and theoreti- 
cal distribution for four different aperture radii. Up- 
per panel: theoretical distribution (red dashed line) and 
observed distributions (solid line histograms). Lower 
panel: difference between theoretical and observed 
distributions. Only objects with detections in 5 epochs 
are included. 
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Fig. 4. — Comparison between theoretical and ob- 
served distributions for a 0.12"aperture in different 
flux bins. Upper panel: theoretical distribution (red 
dashed line) and observed distribution (solid line his- 
tograms). The inset plot show the histogram of catalog 
fluxes, the vertical lines denote the left limit flux of 
the histogram of the given color. Lower panel: dif- 
ference between theoretical and observed distributions 
for same flux bins as in the upper panel. 



Fig. 5. — The C distribution (0.36"aperture) for four 
diff'erent flux bins before (left) and after (right) the 
flux-correction is applied. The solid-line histograms 
show the histograms for the different flux bins, the 
dashed red line shows the theoretical distribution. 
Lower panels shows the differences between the the- 
oretical and observed distributions. For the meaning 
of the histogram colors, see legend and inset plot in 

Fig. a 
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Fig. 6. — Determining the flux-dependent error cor- 
rections for the observed C-distribution. Shown is the 
observed peak of the C-distribution versus the mean 
(red -I- symbols) and median (blue x symbols) flux in 
each bin. The best fit to the data is shown as a dashed 
red line for the mean and a blue dotted line for the me- 
dian. The big plot plot shows the data for low fluxes, 
the small inset shows the same data for the entire flux 
range. The axes labels for the inset plot are identical 
to the big plot. 
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Fig. 8. — Location of the variable (big red stars) and 
non- variable (small black dots) objects in the GOODS 
South field. Objects that have been rejected due to in- 
fluences of nearby bright stars are plotted as small blue 
circles. Only objects with detections in all 5 epochs are 
included. 
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Fig. 7. — Location of the variable (big red stars) and 
non-variable (small black dots) objects in the GOODS 
North field. Objects that have been rejected due to in- 
fluences of nearby bright stars are plotted as small blue 
circles. Only objects with detections in all 5 epochs are 
included. 




Fig. 9. — Sextrator stellarity parameter ClassStar plot- 
ted against z band catalog magnitude. Non-variable 
objects are shown as small black dots, variable objects 
are shown as big red stars. Objects that have been re- 
jected due to the influence of nearby bright stars are 
not shown. 
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Fig. 10. — Flux corrected variability estimator C plot- 
ted against object z catalog magnitude. The color of 
the data point show the Sextractor parameter ClassStar 
which gives an estimate of the stellarity of the object. 
Objects that have been rejected due to the influence of 
nearby bright stars are not shown. 
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Fig. 12. — Variability Strength V in per cent over cat- 
alog magnitude for the final sample of 156 variable 
sources. 
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Fig. 1 1 . — Raw (solid line) and false positive corrected 
(dashed line) number of detections for different mag- 
nitude bins and significance levels. 
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Results of statistical simulations 



Table 4 

Detection power (in per cent) for different mock light curves at a significance of 99.9 per cent with 'erroneous' 
ERRORS. Variability strength V is identical to the one used in Table|2]and defined in the text. Results are given 

FOR an error in THE ERROR MEASUREMENT OF 5 AND 25 PER CENT. . 





Signalstrength 


Slope 5% 


Slope 25% 


Burst 5% 


Burst 25% 


Noise 5%. 


Noise 25%) 


F 


0.5 


0.13 


0.22 


0.15 


0.25 


0.00 


0.00 


C 


0.5 


0.14 


0.26 


0.15 


0.34 


0.00 


0.00 




0.5 


0.14 


1.05 


0.17 


1.18 


0.00 


0.01 


F 


LO 


0.20 


0.34 


0.37 


0.56 


0.11 


0.18 


C 


1.0 


0.20 


0.43 


0.37 


0.69 


0.11 


0.23 


)? 


1.0 


0.22 


1.32 


0.40 


1.92 


0.12 


0.97 


F 


1.5 


0.40 


0.59 


0.98 


1.21 


8.72 


8.75 


C 


1.5 


0.41 


0.76 


0.99 


1.47 


8.80 


9.68 




1.5 


0.45 


2.03 


1.04 


3.49 


8.94 


13.52 


F 


2.0 


0.80 


0.98 


2.43 


2.92 


32.92 


32.17 


C 


2.0 


0.81 


1.20 


2.47 


3.43 


33.07 


33.80 




2.0 


0.87 


2.98 


2.61 


6.50 


33.31 


38.29 


F 


2.5 


1.60 


1.95 


5.71 


6.29 


56.23 


54.96 


C 


2.5 


1.62 


2.33 


5.79 


7.23 


56.36 


56.57 




2.5 


1.70 


4.76 


5.99 


11.62 


56.50 


59.90 


F 


3.0 


3.05 


3.55 


12.16 


12.57 


72.62 


71.01 


C 


3.0 


3.10 


4.17 


12.29 


14.14 


72.73 


72.22 


x" 


3.0 


3.20 


7.47 


12.57 


19.65 


72.82 


74.52 


F 


3.5 


5.56 


6.19 


22.90 


22.13 


82.45 


81.42 


C 


3.5 


5.63 


7.08 


23.10 


24.23 


82.53 


82.31 


x" 


3.5 


5.77 


11.36 


23.34 


30.04 


82.59 


83.83 


F 


4.0 


9.50 


9.74 


36.92 


35.31 


88.55 


87.62 


C 


4.0 


9.62 


11.11 


37.19 


37.92 


88.59 


88.24 


x" 


4.0 


9.88 


16.40 


37.45 


43.16 


88.65 


89.31 


F 


4.5 


15.58 


15.65 


53.77 


50.89 


92.46 


91.50 


C 


4.5 


15.73 


17.30 


54.06 


53.67 


92.50 


91.94 


x" 


4.5 


16.00 


23.46 


54.14 


57.53 


92.51 


92.64 


F 


5.0 


23.68 


23.12 


69.73 


66.09 


94.58 


94.18 


C 


5.0 


23.86 


25.31 


69.96 


68.56 


94.60 


94.50 


x" 


5.0 


24.18 


31.94 


69.95 


70.40 


94.66 


94.99 


F 


5.5 


33.49 


31.93 


82.58 


79.40 


96.24 


95.81 


C 


5.5 


33.73 


34.54 


82.75 


81.21 


96.26 


96.04 


x" 


5.5 


34.25 


41.10 


82.72 


81.52 


96.29 


96.39 


p 


6.0 


45.00 


42.66 


91.52 


88.97 


97.27 


97.02 


c 


6.0 


45.26 


45.37 


91.61 


90.22 


97.28 


97.19 


x" 


6.0 


45.61 


51.54 


91.48 


89.66 


97.28 


97.44 


F 


6.5 


56.66 


53.90 


96.50 


94.88 


97.96 


97.79 


C 


6.5 


56.93 


56.63 


96.54 


95.57 


97.97 


97.90 


x" 


6.5 


57.28 


61.89 


96.48 


94.89 


97.96 


98.07 


F 


7.0 


68.09 


64.65 


98.72 


97.92 


98.39 


98.33 


C 


7.0 


68.33 


67.20 


98.74 


98.21 


98.40 


98.42 


x" 


7.0 


68.55 


71.43 


98.72 


97.72 


98.40 


98.54 


F 


7.5 


78.02 


74.83 


99.66 


99.27 


98.77 


98.69 


C 


7.5 


78.20 


76.94 


99.67 


99.39 


98.77 


98.78 


x" 


7.5 


78.31 


79.91 


99.64 


99.10 


98.77 


98.90 


F 


8.0 


85.96 


82.80 


99.92 


99.78 


99.08 


98.94 


C 


8.0 


86.11 


84.43 


99.93 


99.82 


99.08 


99.01 


X^ 


8.0 


86.16 


86.35 


99.92 


99.67 


99.08 


99.09 


F 


8.5 


91.73 


89.12 


99.98 


99.95 


99.25 


99.18 


C 


8.5 


91.82 


90.28 


99.98 


99.96 


99.26 


99.22 


x" 


8.5 


91.81 


91.57 


99.98 


99.91 


99.25 


99.31 


F 


9.0 


95.46 


93.57 


100.00 


99.99 


99.40 


99.35 


C 


9.0 


95.52 


94.37 


100.00 


99.99 


99.40 


99.38 


x" 


9.0 


95.52 


94.99 


100.00 


99.97 


99.41 


99.45 


F 


9.5 


97.70 


96.39 


100.00 


100.00 


99.48 


99.44 


C 


9.5 


97.72 


96.87 


100.00 


100.00 


99.48 


99.47 




9.5 


97.75 


97.22 


100.00 


100.00 


99.49 


99.52 
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Variable object catalog 



Table 5 

First part of variable objects catalog. Columns: ID: lAU object identification; RA (J2000): Right Ascension 
(Julian Year 2000); Dec (J2000): Declination (Julian Year 2000); mag (z): Sextractor z band magnitude; CS: 

STELLARITY (1 IS POINT LIKE, IS EXTENDED); C: FLUX CORRECTED C FOR 0.36 "aPERTURE; V (%): VARIABILITY STRENGTH IN 

percent; Catalog: Normal (99.9% insignificance) or Clean (99.99% significance); Comment: morphology or 
classification, references for classifications are given; z'- spectroscopic redshift if available and flag (secure (s), 

LIKELY (L), TENTATIVE (T)), 'NoNE' IF SPECTROSCOPY IS AVAIL ABLE BUT NO REDSHIFT WAS FOUND. UNLESS OTHERWISE NOTED, 
ALL SPECTROSCOPIC REDSHIFTS ARE FROm IPoPESSO ET AL.I (120091) ; Mabs'- ABSOLUTE MAGNITUDE IN OBSERVED FRAME Z (ONLY FOR 



OBJECTS WITH SPECTROSCOPI C REDSHIFT > 0). OBJ ECTS THAT HAVE BEEN DETECTED IN OT HER S TUDIES ARE MARKED BY 

SUPERSCRIPTS IN THE ID COLUMN: COHEN ET AL. I (l2006h : C, IKlesman & SARAJEDDd (120071): KS.ITrevese et al.I (120081): T. 



ID 


RA (J2000) 


Dec (J2000) 


mag (z) 


CS 


C 


V (%) 


Catalog 


Comment 




Mabs 


TO'^^^fl^ 00-9749 n f, 


53 01248 


-27 7037879 


25.14 


0.09 


2.55 


28.23 


Clean 


KGSOIVCQ 






J03 3203 .0 1 -274544.7 


53.0125331 


-27.7624232 


24.90 


0.02 


2.05 


22.05 


Normal 


RcsolvGcl 






J033203. 26-274530.3 


53.013574 


-27.7584257 


22.70 


0.03 


2.15 


3.99 


Normal 


Resolved 






J033204.41-274635.5 


53.018367 


-27.7765167 


21.14 


0.98 


3.92 


2.20 




Unresol ved 






J033205. 11-2743 17.5 


53.021302 


-27.721541 1 


21.91 


0.98 


4.30 


3.72 


Clean 


l_Jru"esol ved 






J033205 .40-274429. 2 


53.0224983 


-27.7414518 


22.97 


0.03 


2.47 


7.86 


Clean 


Disk G3.l3.xy 


None 




T03'^90S fiS-974S0S 


S3 03fil S4 


-97 7S9?31 7 


22.70 


0.03 


1.98 


5.01 


Normal 


Resolved 


1 .2964^ 


-22.1 1 


J033209. 57-274634. 9 


53.0398585 


-27.7763722 


22.97 


0.03 


2.03 


4.04 


Normal 


Elliptics.! 






J033209. 58-27424 1.8 


53.0399142 


-27.71 16218 


23.39 


0.03 


2.05 


5.38 


Normal 


Ellipticsl 






J033209. 80-274308. 6 


53.0408277 


-27.7190613 


23.12 


0.03 


1.93 


5.85 


Normal 


Iiitersctin^ system 


2.3021'^ 


-23.20 


J0332 10.52-274628.9 


53.0438334 


-27.7747007 


23.85 


0.03 


2.10 


9.06 


Normal 


Resolved 






J0332 10.9 1-274414.9^^^ 


53.0454704 


-27.7374846 


22.37 


0.93 


4.69 


6.02 


Clean 


Elll|?tlC3l 


1.6082^ 


-23.01 


J03321 1.02-274919.8 


53.045918 


-27.8221721 


23.45 


0.03 


2.15 


9.70 


^Jormal 


InterQctin*^ system 


1.9431-'- 


-23.42 


J033212. 16-274408.8 


53.0506774 


-27.7357915 


24.79 


0.03 


2.29 


22.17 


Clean 


Resolved 






J0332 13.2 1-2747 15.7 


53.0550544 


-27.7876915 


24.09 


0.03 


2.28 


13.89 


Clean 


Resolved 






J033213. 34-274210.5 


53.0555659 


-27.7029202 


23.84 


0.99 


2.34 


7.71 


Clean 


Unresol ved 






J033215. 16-274754.6 


53.0631513 


-27.7985117 


22.98 


0.98 


2.07 


3.32 


Normal 


Unresolved 






J0332 15.93-275329.3 


53.0663591 


-27.8914797 


19.86 


0.92 


2.15 


0.69 


Normal 


Star^' 


0.0000'^ 




J0332 16. 34-27485 1.7 


53 0681031 


-27.8143745 


21.98 


0.98 


2.65 


2.37 


Clean 


Unre solved 






J0332 16.87-2749 16.7 


53.0702817 


-27.8212937 


23.40 


0.98 


2.43 


5.88 


Clean 


Unre solved 


0.0000^ 




J0332 1 7.06-27492 1.9 


53.0710682 


-27.8227401 


19.20 


0.03 


2.03 


0.93 


I'^ormal 


InterQctm*^ system 


033755 


-22.06 


T033917 14.974^03 "X^^-T 


S3 071439fi 


-97 717S8fi4 


20 57 


0.03 


3.14 


1.44 


Clean 


Elliptica.1 






J033217.72-274703.0 


53.0738469 


-27.7841607 


23.59 


0.03 


1.99 


6.18 


Normal 


Interactmg system 






J03321 8.24-27524 1 .4*^^ 


53.0760024 


-27.8781606 


24.28 


0.99 


2.15 


9.61 


Normal 


Unresolved 


2.8049'^ 


-22.57 


J033218.70-275 149.3 


53.0778965 


-27.8637054 


21.74 


0.03 


2.20 


3.04 


Clean 


Disk g3.l3.xy 


0.4568'^ 


-20.55 


J033218.81-274908.5 


53.0783542 


-27.8190409 


23.90 


0.03 


2.55 


14.44 


Clean 


C omple X 






J033218. 84-274529. 2 


53.078519 


-27.7581 146 


18.78 


0.03 


2.44 


0.97 


Clean 


ElliptlC3l 






J033219. 81-275300.9 


53.0825539 


-27.8835727 


24.46 


0.03 


1.95 


11.16 


Normal 


Resolved 


3.7072'^ 


-23.12 


J033219.86-274110.0 


53.0827492 


-27.686119 


23.37 


0.00 


2.07 


8.80 


Normal 


Interacting system 






J033220.80-275144.5 


53.0866794 


-27.8623513 


21.22 


0.98 


2.91 


1.71 


Clean 


Star 






J03322 1.52-274358.7 


53.0896509 


-27.732984 


21.00 


0.97 


1.95 


0.99 


Normal 


StarJ 






J033222.82-274518.4 


53.0950956 


-27.7550986 


23.22 


0.03 


1.92 


5.70 


Normal 


Resolved 






J033223.53-274707.5 


53.0980434 


-27.785425 


23.23 


0.03 


2.00 


7.25 


Normal 


Interacting system 






J033224.23-274129.5 


53.1009433 


-27.691518 


25.01 


0.00 


1.96 


21.51 


Normal 


Complex 






J033224.54-274010.4 


53.1022685 


-27.6695645 


21.88 


0.03 


1.95 


2.54 


Normal 


Interacting system 






J033224.80-274617.9 


53.1033449 


-27.7716431 


23.06 


0.03 


2.69 


7.15 


Clean 


Resolved 






J033225. 10-274403. 2 


53.1045636 


-27.7342138 


17.75 


0.03 


2.73 


0.75 


Clean 


Elliptical 






J033225. 15-274933. 3 


53.1048118 


-27.8259053 


21.61 


0.99 


2.56 


1.73 


Clean 


Stae-' 


0.0000^ 




J033225.99-274142.9 


53.1082952 


-27.6952603 


24.54 


0.03 


1.95 


14.06 


Normal 


Resolved 


0.0459^ 


-12.18 


J033226.40-275532.4 


53.109991 


-27.9256529 


24.32 


0.00 


2.21 


25.29 


Clean 


Complex 






J033226.49-274035.5'^^-'' 


53.1103938 


-27.6765399 


19.60 


0.93 


13.45 


4.35 


Clean 


Core + spiral structure 


0.5404^ 


-22.87 


1033227.01-274105.0*^^ 


53.1125287 


-27.6847238 


19.00 


0.93 


3.97 


1.00 


Clean 


Core + extended emission 


0.7423^ 


-24.31 


J033227.18-274416.5 


53.113269 


-27.73791 


19.37 


0.03 


2.26 


0.89 


Clean 


Elliptical 






J033227.5 1-2756 12.4 


53.114644 


-27.9367764 


24.04 


0.03 


4.18 


20.42 


Clean 


Resolved 


0.6630^ 


-18.97 


J033227.86-275335.6 


53.1160832 


-27.8932186 


21.84 


0.98 


5.00 


2.44 


Clean 


Star' 






J033227.87-275335.9 


53.1161233 


-27.8933035 


21.63 


0.99 


2.30 


1.08 


Clean 


Star' 






J033228.30-274403.6 


53.1179209 


-27.7343234 


23.76 


0.00 


2.25 


28.74 


Clean 


Complex 






J033228.45-274203.8 


53.118527 


-27.7010451 


24.23 


0.99 


2.08 


9.44 


Normal 


Unresolved 






1033229.88-274424.4*^^ 


53.1244949 


-27.7401248 


16.45 


0.03 


2.23 


1.09 


Clean 


Disk galaxy 






1033229.98-274529.9*^^ 


53.1249148 


-27.7583013 


21.06 


0.93 


4.38 


2.62 


Clean 


Core + extended emission 


1.2091^ 


-23.55 


J033229.99-274404.8 


53.1249588 


-27.7346753 


16.84 


0.03 


2.32 


0.84 


Clean 


Disk Galaxy 


0.0746^ 


-20.81 


1033230.06-274523.5*^* 


53.1252547 


-27.756535 


21.81 


0.25 


2.47 


3.20 


Clean 


Resolved 






1033230.22-274504.6*^* 


53.1258995 


-27.7512749 


21.50 


0.92 


5.39 


4.07 


Clean 


Core + extended emission 






J033230.36-275133.2 


53.1264805 


-27.8592312 


25.01 


0.03 


2.09 


22.48 


Normal 


Resolved 






J033231.80-275110.4 


53.1324925 


-27.8528853 


21.03 


0.98 


6.58 


2.99 


Clean 


Unresolved 






J033231.82-275110.6 


53.1325905 


-27.8529435 


21.06 


0.99 


5.74 


2.50 


Clean 


StarJ 






J033231.94-274531.3 


53.1330882 


-27.7587052 


21.11 


0.99 


2.15 


1.16 


Normal 


Star^ 






J033232.04-274523.9 


53.1334996 


-27.7566329 


23.45 


0.03 


1.97 


6.95 


Normal 


Interacting system 






J033232. 12-275636.8 


53.1338346 


-27.9435593 


21.04 


0.99 


1.97 


1.02 


Normal 


Star" 
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Table 5 — Continued 



ID 


RA (J2000) 


Dec (J2000) 


mag(z) 


cs 


C 




Catalog 


Comment 


z 




JUJiziz. Jz-Z /4j Id.4 


53.1346/ /o 


7-7 -77 1 7 1 OA 

-z/. /zlzl89 


22.11 


0.98 


2.24 


2.16 


Clean 


Unresolved 






JUjjZjZ.4y-Z / JU44.U 


C^ 1 7C^/^S7 
7J. i7J70o / 


77 C/155/171 
-Z /.04jD47i 


77 /17 
ZJ.4Z 


U.U7 


7 1 7 
Z. IZ 


O.lj 


Normal 


Disk galaxy 




■■■ 




57 175QQ7^ 

jj.ljjSozo 


77 QQ70^7/; 

-z/.oo/yo3o 


22.93 


0.03 


2.32 


5.00 


Clean 


Interacting system 


u.ys/j 


-2,1.14 






77 fl7Qn4S 
-Z/.0ZyU46 


1 Q 7Q 

10. zy 




7 Q5 

z.yj 


\J.OJ 


Clean 


Elliptical 






JUjJZJ:>.Oo-Z /4UJJ.0 


57 l/in77Q7 


77 ^i7^i54CQ 
-Z/.O /Oj4oy 


7/) Q4 
Z4.y4 


n n7 
u.uz 


7 5Q 

z.Dy 




Clean 


Resolved 






Tn^^O'l'^ ^(i 974704 ^ 




77 7S451 55 


25 21 


U.U4 


1 07 

i.y / 


ZU.Do 


Normal 


Resolved 






juj jZju.yz-z / jjuo.4 


57 1 578 1 07 


77 ftfi5^i70 


24 41 


03 


1 05 

i.yj 


1 1 75 


Normal 


Resolved 






Tn'^^9^7 1^^ 77417k 7 

JUjjZj/.IUZ/HIZO.Z 




77 ^01 1 7^S 


23 36 


03 


2 49 


7 93 


Clean 


Resolved 






Tnn7T7 77/iA;nQ iC 
JUJ7ZJ / .yj-z /4Duy. i 


57 1 5ttA7/15 


77 KM 1 07 

-z /. /oy lyz 


ly.yo 


n 07 


A 77 
4.7 / 


J .J 1 


Clean 


Unresolved 


U.UoO 


-lo.UZ 


JUJ7Z75. IZ-Z /7y44.5 


57 1 5QQ7n7 
J7. 1 J607U/ 


77 A^7/1/l/1/l 
-Z / .ODZ4444 


7n A'X 


A 07 

u.yz 


1 7 AO 

1 j.uy 


5.39 


Clean 


Core + extemded emission 


\J.OJ IK) 


01 on 

-Zj.ZU 


JU77Z70.oy-Z / J4UD. / 


57 1 /;7rt5/15 


77 Qni Q5ni 
-z/.yuiojul 


7/1 /IQ 

Z4.4y 


n n7 

U.U7 


1 OO 

i.yy 


18.00 


Normal 


Resolved 






Tm777Q nQ 77y)/;ni qCSTS 


57 1/;7Q5Q7 


77 7<;7 1 fJV^ 

-2.1 .1131 


20.96 


0.94 


7.88 


4.39 


Clean 


Core + extended emission 






TnQ777D /17 77?^nn 5 


57 1 /^y) /I 5/^7 
j3. 104430/ 


77 QQ7/1/^QQ 

-2/.iSo34ooy 


20.54 


0.03 


2.01 


1.23 


Normal 


Elliptical 


"■ 




J [JiiJAV.Z l-L /4y4y. / 


57 1 /;7Qm 7 
jJ.iO/oUl / 


77 Q7n/tQ1 

-z/.is3U4ol 


24.69 


0.02 


1.99 


17.90 


Normal 


Complex 


None 




JU37z4U.5y-z /j44y.z 


57 1 7n7/;7Q 


77 Q1 7/^/^/17 

-z/.yi30043 


25.14 


0.02 


2.01 


17.08 


Normal 


Resolved 






Tm07/11 AC 7'7C71/1 1 
JUJ3Z41.U7~Z/7Z74.1 


51 171 AC Al 

53.1 /1U54/ 


7-7 CViC 1 /I 

-z/. 8/61468 


20.72 


0.97 


2.71 


1.08 


Clean 


Star^ 


u.uUuU 




J03324 1 .87-27465 1 . 1 


53.1744478 


-27.7808655 


23.35 


0.01 


2.01 


12.07 


Normal 


Tadpole galaxy"^ 






J033242.6 1 -275453.8 


53.1775348 


-27.9149453 


20.82 


0.99 


2.64 


1.25 


Clean 


Star^ 


0.0000* 




J()33243. 24-2749 14. 2'''' 


53.1801493 


-27.8206046 


22.49 


0.93 


9.56 


12.41 


Clean 


Core + extended emission 


0.2145^ 


-17.66 


J033243. 93-27435 1.1 


53.1830401 


-27.7308524 


24.26 


0.04 


2.00 


9.75 


Normal 


Resolved 






J033244. 10-275212.9 


53.1837535 


-27.8702564 


20.88 


0.94 


2.33 


0.90 


Clean 


StarJ 






J033245.02-275207.7 


53.1875887 


-27.8688155 


22.57 


0.03 


4.45 


7.33 


Clean 


Disk galaxy 






J033246.37-274912.8 


53.1932061 


-27.8202112 


21.94 


0.03 


2.08 


2.07 


Normal 


Elliptical 


0.6830^ 


-21.15 


J033246.39-274820.1 


53.1932909 


-27.8055737 


21.16 


0.97 


2.59 


1.45 


Clean 


Star^ 






J033247 .53-275 159.9 


53.1980298 


-27.8666421 


21.01 


0.99 


2.11 


0.95 


Normal 


Star 






J033247.98-274855.7 


53.1999289 


-27.8154702 


20.56 


0.04 


2.61 


1.33 


Clean 


Elliptical 


0.2340^ 


-19.79 


J033251.22-275418.3 


53.2134115 


-27.905076 


24.25 


0.03 


1.95 


13.30 


Normal 


Resolved 






J033252.88-275119.8'^' 


53.2203537 


-27.8555099 


21.84 


0.12 


2.93 


4.23 


Clean 


Disk galaxy 


1.2283^ 


-22.81 


J033253.44-275001.4 


53.2226606 


-27.8337103 


24.44 


0.00 


1.96 


15.49 


Normal 


Complex 






J123553. 12-1-621037.5 


188.9713311 


62.1770907 


21.26 


0.93 


5.71 


4.78 


Clean 


Interacting system 






JI23556.88+62I1I7.3 


188.9870032 


62.1881435 


20.23 


0.92 


6.41 


2.34 


Clean 


Unresolved 






J123557.62+621024.7 


188.9900912 


62.1735292 


22.62 


0.94 


4.55 


6.53 


Clean 


Core + extended emission 






J123603.82+621039.3 


189.0159143 


62.1775909 


24.48 


0.00 


2.60 


21.88 


Clean 


Interacting system 






J123605.75+621356.1 


189.0239559 


62.2322631 


20.82 


0.99 


4.01 


1.93 


Clean 


Unresolved 






J123606.47+621506.4 


189.0269558 


62.2517698 


21.42 


0.89 


1.98 


1.28 


Normal 


Core + extended emission 






J123612.6I+621238.4 


189.0525622 


62.2106579 


20.18 


0.89 


3.11 


1.05 


Clean 


Core + extended emission 






JI23617.99-f621635.3 


189.0749736 


62.2764767 


20.18 


0.87 


9.25 


4.21 


Clean 


Core + extended emission 






JI23619.57-I-620715.2 


189.0815546 


62.1208771 


25.14 


0.02 


2.04 


22.65 


Normal 


Resolved 






JI23621.26-I-621640.4 


189.0885649 


62.2778893 


21.57 


0.03 


1.94 


1.88 


Normal 


Interacting system 






J123622.94+62I527.0 


189.095583 


62.2574929 


20.27 


0.90 


4.52 


1.70 


Clean 


Core + extended emission 






J123627.48+621406.4 


189.1144932 


62.2351164 


25.22 


0.04 


2.25 


25.05 


Clean 


Resolved 






J123627.98+621508.1 


189.1165918 


62.2522625 


21.86 


0.03 


1.94 


2.33 


Normal 


Interacting system 






J123629.44+62I5I3.3 


189.1226602 


62.2536973 


23.70 


0.98 


1.97 


5.83 


Normal 


Um-esolved 






J123629.68+62I734.7 


189.1236482 


62.2929798 


22.45 


0.98 


2.27 


2.73 


Clean 


Unresolved 






J123631. 68-1-620848.7 


189.1320148 


62.1468517 


22.56 


0.03 


2.14 


3.24 


Normal 


Elliptical 






J123631. 70-1-620752.3 


189.1320712 


62.1311963 


24.12 


0.03 


2.19 


9.17 


Clean 


Resolved 






J123632.50-I-620701.7 


189.1354166 


62.1171419 


25.37 


0.82 


2.03 


20.07 


Normal 


Resolved 






J123633. 23-1-620834.9 


189.1384756 


62.1430369 


21.08 


0.03 


5.24 


3.87 


Clean 


Interacting system 






JI23636.5I4-620806.4 


189.1521369 


62.1350985 


22.29 


0.03 


2.06 


3.24 


Normal 


Elliptical 






J123637.85+620724.1 


189.1576928 


62.1233725 


22.21 


0.00 


1.94 


5.06 


Normal 


Complex 






J123641.44+620730.4 


189.1726778 


62.1251072 


22.62 


0.00 


1.94 


5.58 


Normal 


Disk galaxy 






J123648.32-I-621250.1 


189.2013266 


62.2139299 


20.39 


0.93 


2.93 


1.02 


Clean 


Core + extended emission 






J 123650.44-1-620749.7 


189.2101509 


62.1304753 


23.62 


0.03 


1.98 


6.34 


Normal 


Resolved 






J123650.75-I-621439.9 


189.2114589 


62.2444226 


23.73 


0.03 


2.13 


9.86 


Normal 


Disk galaxy 






J123651.32-I-621751.2 


189.2138382 


62.2975491 


24.46 


0.03 


1.98 


13.91 


Normal 


Resolved 






J123652.44-I-620959.9 


189.2185084 


62.1666296 


22.44 


0.98 


2.23 


2.58 


Clean 


Unresolved 






J123654.99-I-621635.1 


189.2291393 


62.2764217 


25.46 


0.97 


2.31 


22.23 


Clean 


Unresolved 






JI23655.904-620828.3 


189.2329249 


62.1411947 


21.01 


0.99 


2.85 


1.49 


Clean 


Uru"esolved 






J123656.51+620837.7 


189.2354508 


62.1438079 


20.75 


0.99 


3.38 


1.56 


Clean 


Unresolved 
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Table 5 — Continued 



ID 


RA (J2000) 


Dec (J2000) 


mag iz) 


CS 


C 


V{%) 


Catalog 


Comment z ^abs 


J123656. 91 +621950.3 


189.2371093 


62.330639 


21.61 


0.03 


2.07 


1.61 


Normal 


Disk galaxy 


J 123700.7 1 +62 1 854.4 


189.2529566 


62.3151 155 


23.58 


0.03 


2.20 


9.72 


Clean 


Interacting system 


J123700.88+621 129.5 


189.2536832 


62.1915322 


20.32 


0.97 


3.40 


1.22 


Clean 


Star'' 


J123701.55+622103.9 


189.2564508 


62.3510871 


23.02 


0.06 


2.04 


3.85 


Normal 


Elliptical 


J123702.09+621737.8 


189.2586984 


62.2938331 


20.43 


0.99 


2.02 


0.79 


Normal 


Umesolved 


J123702.72+621543.9 


189.2613473 


62.2621972 


19.41 


0.03 


2.18 


0.91 


Clean 


Disk galaxy 


J123704.80+621455.2 


189.2699968 


62.2486665 


20.99 


0.98 


2.79 


1.44 


Clean 


Umesolved 


J123705.48+621526.8 


189.2728214 


62.257434 


24.62 


0.03 


2.16 


14.31 


Normal 


Resolved 


J123706.25 +622136.9 


189.2760529 


62.360261 1 


20.90 


0.90 


10.05 


8.36 


Clean 


Disk galaxy 


J123706.87+621702.5 


189.2786271 


62.2840158 


19.86 


0.92 


3.81 


1.27 


Clean 


Core + extended emission 1.0200^^ -24.29 


J123706.93+621429.9 


189.2788629 


62.2416441 


20.70 


0.98 


2.77 


1.24 


Clean 


Unresolved 


J123707 .49+622148.1 


189.2812076 


62.3633707 


22.03 


0.92 


3.75 


3.83 


Clean 


Elliptical 


J123708.35+621105.8 


189.284812 


62.184941 


23.47 


0.03 


2.00 


6.04 


Normal 


Resolved 


J123708.65 +621051.5 


189.286038 


62.1809766 


20.67 


0.03 


2.03 


1.91 


Normal 


Interacting system 


J123716.68+621733.6 


189.3194798 


62.2926732 


22.12 


0.72 


2.97 


3.74 


Clean 


Elliptical 


J123717.79+622034.2 


189.3241356 


62.3428338 


23.91 


0.00 


1.94 


26.94 


Normal 


Resolved 


J123717. 82+621 130.5 


189.3242369 


62.1918003 


24.28 


0.17 


2.19 


10.70 


Clean 


Resolved 


J123719.47+621320.4 


189.3311376 


62.2223455 


20.79 


0.96 


2.56 


1.20 


Clean 


Core + extended emission 


J123720.16+621518.9 


189.3339866 


62.2552618 


20.16 


0.88 


6.72 


2.39 


Clean 


Core + extended emission 


J123723.70+621200.3 


189.3487632 


62.2000936 


24.82 


0.00 


1.97 


21.53 


Normal 


Resolved 


J123723.72+622113.3 


189.3488236 


62.3537049 


23.62 


0.98 


2.30 


6.61 


Clean 


Unresolved 


J123724.77+622103.0 


189.3531954 


62.3508303 


23.80 


0.00 


1.96 


14.08 


Normal 


Resolved 


J123728.43+622044.8 


189.3684533 


62.3457887 


22.69 


0.02 


2.24 


7.50 


Clean 


Resolved 


J123728.95+621 127.8 


189.3706346 


62.1910553 


21.75 


0.03 


6.06 


6.20 


Clean 


Supernova'^ 


J123729.58+621557.8 


189.3732355 


62.2660663 


19.91 


0.85 


2.12 


0.60 


Normal 


Core + extended emission 


J123732.41+621751.4 


189.3850213 


62.2976196 


24.58 


0.03 


2.55 


17.65 


Clean 


Resolved 


J123736.11+621619.1 


189.4004397 


62.2719614 


24.07 


0.03 


2.91 


15.16 


Clean 


Resolved 


J123736.59+621632.7 


189.4024445 


62.2757619 


25.03 


0.03 


3.09 


37.46 


Clean 


Resolved 


J123738. 83+622024.0 


189.4117843 


62.3399951 


20.20 


0.96 


3.41 


1.24 


Clean 


Core + extended emission 


J123740.64+622007.9 


189.4193208 


62.3355268 


22.60 


0.03 


2.50 


5.41 


Clean 


Disk galaxy 


J123741.21+621925.2 


189.4217126 


62.3236767 


21.04 


0.98 


2.84 


1.50 


Clean 


Unresolved 


J123741.38+621540.2 


189.4224198 


62.2611689 


25.19 


0.02 


2.11 


20.13 


Normal 


Resolved 


J123742.12+621903.0 


189.425516 


62.3174941 


24.48 


0.04 


3.43 


29.44 


Clean 


Resolved 


J123742.53+621812.2 


189.4272043 


62.3033909 


21.31 


0.98 


5.40 


3.36 


Clean 


Unresolved 


J123746.85+621624.2 


189.4452252 


62.2733875 


21.94 


0.98 


2.27 


1.96 


Clean 


Unresolved 


J123749.58+621346.6 


189.4565982 


62.2296055 


24.56 


0.98 


2.49 


14.73 


Clean 


Unresolved 


J123754.25+621853.0 


189.47606 


62.3147251 


24.84 


0.02 


2.21 


20.84 


Clean 


Resolved 
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